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1. Introduction. 

This note presents an attempt to provide a conceptual framework for variational formulations of 
classical physics. It is based on a study of static systems in collaboration with P. Urbanski of Warsaw. 

Variational principles of physics have all a common source in the principle of virtual work well known 
in statics of machanical systems. This principle is presented here as the first step in characterizing 
local stable equilibria of static systems. An extended analysis of local equilibria is given for systems 
with configuration manifolds of finite dimensions. Numerous examples of the principle of virtual work 
and the Legendre transformation applied to static mechanical systems are provided. Configuration 
spaces for the dynamics of autonomous mechanical systems and for statics of continua are constructed 
in the final sections. These configuration spaces are not differential manifolds. 

The final section of the present note is an initial part of a programme of geometric formulations of 
the foundations of statics of continua undertaken jointly with F. Cardin of Padova in collaboration 
with P. Urbanski. 

2. Configurations, processes, work. 

Experimental study of physical systems consists in subjecting a system to external control and 
observing the results. A theoretical description of the system should have the power to predict observed 
results. The formulation of a theoretical description must start with assigning to the system parameters 
characterizing its configuration. These parameters do not ever exhaust all characteristics of the system. 
The choice is dictated by the interactions to which the system will be subjected. Statics of mechanical 
systems assigns to a mechanical system points in a differential manifold characterizing the geometrical 
placement of the system in the physical space. The points are translated in parameters by means of 
coordinate systems. The differential manifold is called the configuration space. Control of a static 
mechanical system consists in coupling the system with external devices with the same configuration 
space and observing the resulting configurations. At least two configuration spaces are involved. The 
control configuration space represents degrees of freedom which are being controlled. The observed 
configurations are elements of the observed configuration space. The two spaces are not necessarily the 
same. The configuration spaces considered represent configurations of stable local equilibrium. 

It is sometimes necessary to redefine the control configuration space. It is sometimes found that the 
same control arrangements do not result in the same response. When this happens we must discover 
and include in the description additional hidden degrees of freedom participating in the process of 
control. The configuration space with these degrees of freedom included is the internal configuration 
space. We will assume initially that the inclusion of additional degrees of freedom is not necessary and 
that the internal configuration space is the same as the control configuration space. 

An important form of control consists in subjecting a system to a quasi static process by coupling it 
to an external device with slowly changing configuration. The changes must be slow enough to justify 
the treatment of a process as a succession of configurations of static equilibrium configurations rather 
than a dynamic event. Criteria of equilibrium are based on the response to quasi static processes. 
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Quasistatic processes are represented as oriented embedded arcs in the configuration space. Let 
Q be a differential manifold of dimension m representing the configuration space of a static system. 
An embedded arc is a subset c C Q which is the image of an embedding 7: [0, a] — + Q of a closed 
interval [0,a] C K. defined as the restriction to [0, a] of an embedding 7:R — > Q. The embedding 
7 is a representative of an embedded arc called its parameterization. Different embeddings may be 
parameterizations of the same arc. A parameterization 7: [0, a] — > Q of an arc c induces an orientation 
of the arc. The arc is oriented from 7(0) to 7(a) if the two points are distinct. If 7(0) = 7(a), then 
the arc is a cycle. We will exclude cyclic processes from the discussion. They are of no interest for 
developing criteria of equilibrium. Two parameterizations 7 and 7' induce the same orientation if 
7'(0) = 7(0). The boundary dc of an arc c is a set of two points q and q' such that 7(0) = q and 
7(a) = q' for some parameterization 7: [0, a] — > Q of c. The designation of one of the boundary points 
as the initial configuration of the process represented by an arc specifies an orientation. An oriented 
embedded arc is defined as a pair (c, go), where c is an embedded arc and qo is a point in the boundary 
dc designated as the initial configuration. The remaining boundary point is the terminal configuration 
of (c, qo). A parameterization 7: [0, a] — > Q of an arc c such that 7(0) = qo is said to be compatible 
with the orientation of the oriented arc (c, qo). Compatible parameterizations will be used for oriented 
arcs. 

Let (c, q ) be an oriented arc. We refer to points in c as elements of (c, q ). Elements of an oriented 
arc (c, q ) are ordered by the relation ^ defined in terms of any parameterization 7 compatible with 
the orientation. The relation j(s') 7(3) is equivalent to s' ^ s. Relations <, ^, and > are defined 
in a similar way. They reflect the corresponding relations between the values of the parameter. 

Let (c, qo) and (c', q' ) be oriented arcs such that c' C c. There are two ordering relations in (c', q' ) 
since elements of (c',g ) are also elements of (c, qo). The arc (c',q' Q ) is said to be included in the arc 
(c, g ) if the two ordering relations coincide. The inclusion relation is denoted by (c', q' ) C (c, qo). The 
process represented by {c',q' Q ) is a subprocess of the process represented by (c, go) if ( c '>9o) c ( c :<Zo)- 
If (c 1 , <7o) an d (c 2 , g 2 ) are oriented arcs and if one of the pairs (c 1 Uc 2 , gg) or (c 1 Uc 2 , gg) is an oriented 
arc, then it is considered the union of (c 1 , gg) and (c 2 , g 2 ) and is denoted by (c 1 , gg)U(c 2 , g 2 ). Similarly 
if one of the pairs (c 1 n c 2 , gg) or (c 1 n c 2 , g 2 ) is an oriented arc, then it is considered the intersection 
of (c^gg) and (c 2 ,g 2 ) and is denoted by (c^gg) n (c 2 ,gg). 

Let (c, go) be an oriented arc. A function a: c — > K is considered a function on (c, go). This function 
is said to be differentiable if it the restriction to c of a differcntiablc function on Q. Let AQ be the 
algebra of differentiable functions on Q and let lo(Q, c) C AQ be the ideal of functions vanishing on c. 
Two functions / and /' on Q have the same restrictions to c if and only if /' — / € \o{M, S). It follows 
that the algebra A(c, q ) of differentiable functions on (c, q ) is canonically isomorphic to the quotient 
algebra AQ/\ (Q, c). The orientation of (c, g ) makes it possible to single out increasing functions. A 
function h G Ac is increasing if q' > q implies h(q') > h(q). 

An oriented arc represents a quasi static process. We use these two terms interchangeably. We 
denote by PQ the set of all processes in Q. An essential part of the characterization of a static system 
is the specification of a set C C PQ of admissible processes in the configuration space Q. Admissible 
processes are the processes that can be actually induced using the control devices at our disposal. The 
following conditions are satisfied. 

(1) A subprocess of an admissible process is admissible. 

(2) The union of admissible processes is admissible. 

(3) If a process can be approximated with admissible processes, then it is admissible. 

The last condition is made more precise by the following statement. Let g" l7 g" 2 , . . . , g" i; . .. be a 
sequence of interior configurations in a process (c, go) converging to the terminal configuration gi. 
If processes (c q t.\,qo) are admissible, then (c, go) is admissible. The symbol (c 9 |,g ) denotes the 
subprocess of (c, g ) with g e c as the terminal configuration. It is based on the arc 

c q \ = {q 1 £c; q' ^q}. (1) 
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Another object characterizing a static system is a function 



W:C^R (2) 

associating with each admissible process (c, q ) the work W(c, q ) of the process. Since subprocesses of 
admissible processes are admissible we can associate with each admissible process (c, go) the function 

f W(c g \,q ) if q + qo 
w { ^ qo) :c^R:q^ i (3) 
10 if q = q . 

The work function will be assumed to satisfy the following conditions. 

(1) Work is additive in the sense that if (c, q Q ) is the union of admissible processes 

(c%8),(cS ft 1 ),. ■■,(<:",«?) (4) 

such that for i = 0, 1, . . . , n — 1 the terminal configuration of (c l , g ) is the initial configuration 
of (c i+1 ,g +1 ), then 

71 

W(c,q )=J2W(c i ,ql). (5) 

i=0 

(2) For each admissible process (c, go) the function iO( c ,g ) is differentiable. 

3. An alternative definition of orientation. 

For the purpose of constructing an orientation of an embedded arc c we introduce the tangent 
bundle Tc. Let 7: [0, a] — > Q be a parameterization of c and let 7: M — > Q be an embedding such 
that 7 = 7|[0, a]. The image c = 7(E) is a submanifold of Q and has a well defined tangent bundle 
Tc C TQ. The tangent bundle Tc is the restriction 

Tc = TcnrQ 1 (c) (6) 

of TC to c. The result of this construction is independent of the choice of the parameterization. An 
orientation of the arc c is a subset o € Tc with the following properties. 

(!) T Q(°) = c - 

(2) If v G o then kv e o for each k > 0. 

(3) If w £ o then kv ^ o for /c 0. 

(4) There is a parameterization 7: [0, a] — > Q such that t7(s) £ o for each s £]0, a[. 

The symbol t7(s) denotes the vector tangent to the curve 7 at 7(5). There are obviously two different 
orientations of each embedded arc. An oriented embedded arc is a pair (c, o) of an embedded arc c 
and one of its orientations o. Parameterizations of an arc satisfying the condition (4) above are said 
to be compatible with the orientation o. A parameterization 7: [0, a] — > Q can be used to generate an 
orientation of an arc. The set 

o = {v G Tc; 3 se[0;O]ife>0 « = fct7(*)} (7) 

is an orientation of the arc c = 7([0, a]) and the parameterization 7 is compatible with this orientation. 
Vectors t7(0) an Vy(a) are obtained as limits. Specifying the initial configuration go determines an 
orientation o of the arc unless the terminal configuration is coincides with the initial configuration. If 
such are excluded, then the pair (c, go) is a representation of the oriented arc equivalent to (c, go). 
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4. Stable local equilibrium configurations. 

Let C° C Q be the set of initial configurations of all admissible processes for a static system. A point 
g G C° is called a stable local equilibrium configuration of the system if for each admissible process 
(c, go) initiating at go the function u>( c ,<j ) nas a l° ca l minimum at g . 

The function w>( Ci9o ) will have a local minimum at go if there is a point gee such that the function 
w (c.q ) restricted to the arc 

C q \ = {q' ec; q' ^ q} (8) 

is increasing. This observation implies that equilibrium at go is a property of the germ of W( Ci9o ) at go. 
We will state the definition of equilibrium in terms of germs. 

In the set PQ of oriented embedded arcs in Q we introduce an equivalence relation. Arcs (c, g ) 
and (c',g ) are equivalent if g = go and if there is an open neighbourhood U C Q of go such that 
c' n U = c n U. Equivalence classes are called germs of processes. The set of germs is denoted by P C Q. 
The germ of a process (c, g ) is denoted by j c (c, g ). Germs of processes can be defined in terms of 
ideals. In the algebra AQ we introduce an ideal l c (Q,g) for each g G Q. A function / is in \ c (Q,q) if 
there is a neighbourhood U C Q of q such that f\U — 0. Arcs (c, go) and (c', q' Q ) are equivalent if 

\c(Q, q' ) + lo(Q, c') = l c (Q, g ) + l (Q, c). (9) 

This is the same equivalence relation as the one defined above. The symbol C c qo will denote the set of 
germs of admissible processes initiating at go- 

With a process (c, go) we associate an ideal l c (c, go) in the algebra A(c, go) of diffcrcntiable functions 
on (c, go). A function h is in this ideal if there is a neighbourhood V G c of the initial configuration go 
such that h\V = 0. The algebra A(c, go) is canonically isomorphic to the quotient algebra AQ/\ (Q, c) 
and the ideal l c (c, go) is isomorphic to the quotient 

(lc(g,go) + lo(0,c))/l (0,c). (10) 

We will adopt the following identification 

A c (c,g ) = A(c,g )/I c (c,g ) is identified with AQ/ (\ C (Q, g ) + \ (Q, c)) . (11) 

As a consequence of this identification the algebra A c (c, g ) is associated with the germ j c (c, q ) rather 
than with the process (c, g ) since the algebra A c (c', q' ) is the same as the algebra A c (c, g ) if (c', q' ) 
and (c, go) are equivalent. 

Each function h G A(c, go) has a germ } c h G A c (c, go). The germ j c W( CiIJ0 ) of the work function W( c go ) 
along an admissible process (c, q ) is the same as the germ j c W( c ' if (c', q' ) and (c, g ) are equivalent. 
An element of A c (c, go) is said to be increasing if it has an increasing representative in A(c, go). 

The criterion of stable local equilibrium is based on response of the system to the initial phase of a 
process. It is therefore correct to adopt the following definition of equilibrium formulated in terms of 
germs. 

A point g G C° is called a stable local equilibrium configuration if for each germ j c (c, go) G C^ the 
germ j c tO( c , 9o ) is increasing. 

We have excluded constant processes represented by constant arcs. In a complete discussion of 
equilibrium configurations isolated admissible configurations and admissible constant processes should 
be considered. Only constant admissible configurations can initiate at an isolated admissible configu- 
ration. Such configurations ar obviously stable local equilibrium configurations. No further discussion 
is necessary. 

5. Differential criteria of local equilibrium. 



4 



5.1. Jets of processes and functions. 

We will denote by K the set N U {oo,c}, where c stands for the cardinality of R. The ordering 
relations sC, <, ^, and > have in IK the usual meaning of inequalities of cardinal numbers. Let q be a 
point in a differential manifold Q. In the algebra AQ of differentiablc functions on Q we introduce a 
sequence of ideals 

\ Q (Q,q), h(Q, <?),..., loo(Q,«), lc(Q,«). (12) 

The ideal lo(<3, q) associated with q is maximal in the sense that it is not a proper subset of any ideal 
except the trivial ideal AQ. It is known that all maximal ideals in AQ are associated with points. 
For k G N, the ideal \k(Q,q) is the power (lo(Q, q)) k+1 of the ideal \o(Q,q)- The ideal loo(Q,g) is the 
intersection f] keN h(Q, q)- The ideal l c (<3, q) is the set of functions each vanishing in a neighbourhood 
of q. This ideal was introduced earlier. Inclusion relations 

\ k (Q,q) C \ k >{Q,q) (13) 

hold for all k! and k in K such that k! < k. 

For each k e K we introduce an equivalence relation in the set PQ of oriented embedded arcs in Q. 
Arcs (c, g ) and (c',q' ) are equivalent if 

h(Q, q'o) + \o(Q, c') = h(Q, Qo) + \o(Q, c). (14) 

Equivalence classes are called k-jets of processes. The c-jets are the germs introduced earlier. The set 
of fc-jets will be denoted by P k Q. The fc-jet of a process (c, g ) will be denoted by j fc (c, q ). Inclusion 
relations imply the existence of projections 

7T k ' k q: P k Q - P k 'Q: j fc (c, go) ~ j fc '(c, q ) (15) 

for k! < k in addition to 

7r kQ :P k Q^Q: ] k (c,q )^qo. (16) 
Corresponding to the sequence (12) of ideals we have the sequence 

A°(Q, q) = AQ/\ (Q, q),..., A°°(Q, q) = AQ/UQ, q), A C (Q, g) = AQ/\ C (Q, q) (17) 

of quotient algebras. 

In the algebra A(c, go) of differentiable functions on the oriented arc (c, go) we have ideals 

l (c,g ), li(c,g ), ... , loo(c,g ), l c (c,g ) (18) 

defined just as the ideals in the sequence (12) and the corresponding quotient algebras 

A°(c,g ) = A(c,g )/I (c,g ), A^c.go) = A(c, g )/li(c, g ), . . . 

A°°(c,g ) = A(c,g )/loo(c,g ), A c (c, g ) = A(c, g )/l c (c, g ). (19) 

For each k € K the ideal U(c, q ) is isomorphic to the quotient 

(ifc(Q,«b) + io(Q,c))/l (Q,c). (20) 

This justifies the following identification 

A fe (c,g ) = A(c,g )/lfe(c,g ) is identified with AQ/ (\ k (Q, g ) + \ (Q, c)) . (21) 

The identification implies that the algebra A k (c,q ) is associated with the A; -jet j fe (c, g ) rather than 
with the process (c, q ) since the algebra A k (c',q' ) is the same as the algebra A fe (c, q ) if (c',g ) and 
(c, go) are in the same jet. 

An element of the quotient algebra A fc (c, g ) = A(c, go)/lfe(c, go) is said to be increasing if it has an 
increasing representative in A(c, go). 
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5.2. Jets of functions on [0, a] C R. 

The algebra A[0, a] of differentiable functions on [0, a] C R, its ideals 

l ([0,o],0), li([0,o],0), ... , loo([0,o],0), l e ([0,o],0). (22) 

and the quotient algebras 

A°([0,a],0) = A[0,a]/l ([0,o],0), A\[0, a], 0) = A[0, o]/li([0, a], 0), . . . 

A oc ([0,a],0)=A[0,a]/l oc ([0,a],0), A c ([0, a], 0) = A[0, o]/l e ([0, a], 0) (23) 

deserve special attention. Elements of an algebra A fc ([0, a], 0) are fc-jets of functions on [0, o]. 

A differentiable function on [0, a] is the restriction to [0, a] of a differentiable function on R. A func- 
tion is in an ideal lfc([0, a], 0) with k < c if and only if it is the restriction to [0, a] of a function in U(R, 0). 
This is not true for k = c. As a consequence the quotient algebra A fe ([0,a],0) = A[0, a]/lfc([0, a], 0) is 
isomorphic to A fe (R, 0) = AR/l fe (R, 0) for k < c but not for k = c. In each algebra A fc ([0, a], 0) there is 
a sequence of ideals 

l fc o ([0,a],0) = lo([0,a],0)/l fc ([0,a],0), l fe i([0, a], 0) = I^O, a], 0)/U([0, a], 0), . . . 

l\([0,a],0) = l fe ([0, a ],0)/l fc ([0,a],0). (24) 

The quotient algebra A fc ([0, a], 0)/l fe fe < ([0, a], 0) with k' ^ fe is isomorphic to A fc '([0,a],0). 

A function g G A[0, a] is said to be increasing if the inequality s' > s implies q(s') > g(s). An element 
of a quotient algebra A c ([0, a],0) = A[0, a]/l c ([0, a], 0) is said to be increasing if it has an increasing 
representative in A[0, a]. A function g G A[0, a] is in l c ([0, a],0) if and only if there is a real number 
5 > such that g(s) = for s < <5. It follows that the germ j c <7(0) of a function g e A[0, a] is increasing 
if there is a number <5 > such that g is increasing in [0, 5}. 

PROPOSITION 1. For k e N a function g E A[0,a] is in l fc ([0, a],0) if and onfy if D*5(0) = for each 
i ^ fc. 

Proof: The derivatives of a function g £ A[0, a] at are well defined and the function can be repre- 
sented by the Taylor formula 

g = e (g) + ei(g)a + ... + e k (g)s k + rs k+ \ (25) 

where 

ei(g) = ^D\ 9 (0), (26) 

s: [0, a] — > R is the canonical injection, and r is a differentiable function on [0, a]. The function s is 
in l ([0, o],0) and the power s k+1 is in l fc ([0, a],0). If D l g(0) = for each i ^ k, then g = rs k+1 is in 
l fc ([0,o],0). 

A function g G l;([0, o],0) is a combination of products goffi ' " ■ 9l °f elements of l ( [0, a],0). The 
derivative is a combination of products of functions with each product containing at least I factors 
in l ([0, a],0). It follows that the derivative Dg of a function g G lj([0, a],0) is in lj_i([0, a],0). If 
g G l fc ([0,a],0), then D°.g = g G l fc ([0,o],0), = D.o G l fc _i([0,o],0), B 2 g G l fc _ 2 ([0, a], 0), 
D k g G l ([0,a],0). Hence, D^(0) = for each i < fc. a 

A function .g G A[0, a] is in l°°([0, a], 0) if and only if D 4 .g(0) = for each i G N. 
It follows from the Proposition 1 that for k < c the jet ] k g(0) G A fc ([0, a],0) is fully represented by 
the sequence 

eo(fiO, e i(5), ■ •• - e k(g) (27) 

of derivatives 

e*(5) = ^D^(O) (28) 
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of its representative g G A[0, a]. If k G N, then the polynomial 

eo(5) + ei(5)s + ... + e fe ( 5 )s fc (29) 

represents the jet j fc g(0) and the product of two jets is represented by the product of the corresponding 
polynomials truncated after the first k + 1 terms. 

For k < c an element of the ideal l fc o([0, a], 0) C A fe ([0, o],0) is represented by the sequence 
eo, ei,... , e/j with eo = 0. This element is said to be positive if the first non zero element in 
the sequence is positive. The element is said to be negative if the first non zero element in the sequence 
is negative. Each element of the ideal l fc o([0, a], 0) is either positive or negative if it is not zero. There 
are obvious relations >, <, ^, and < between elements of l fc ([0, a], 0). 

PROPOSITION 2. Let k < c. If a jet ] k g(0) in the ideal l fe ([0, a], 0) is positive, then the germ] c g(0) is 
increasing. If the jet } k g(0) is negative, then the germ] c g(0) is decreasing. 

Proof: If ei{g) is the first non zero element in the sequence (27) of derivatives, then the function g 
is represented by Taylor formula 

g(s) = e l (g)s l + r(s)s l+1 . (30) 



From 



lim ((I + l)r{s)s + Dr(s)s 2 ) = (31) 



it follows that there is a number 5 > such that 

|(Z + l)r( S ) S + Dr( S ) S 2 | < \l ei (g)\ (32) 

for \s\ < S. If ei{g) > 0, then the function g is increasing in the interval [0, 6] since the derivative 

D 5 = leis 1 - 1 + + l)r(s)s l + Dr(s)s l+1 

= (lei + + l)r(*)a + Dr(s)s 2 ) s 1 ' 1 (33) 

is positive for < s < S. This is a consequence of the Lagrange mean value theorem. It follows that 
the germ ] c g(0) is increasing. It is shown in a similar way that if e/ is negative, then the germ is 
decreasing. B 

It follows from the proposition that if the germ j c g(0) of a function g is increasing, then the jet j fc g(0) 
is non negative for each k < c. 



5.3. Conditions of equilibrium. 

Let h be a function in the algebra A(c, qo) and let 7: [0, a] — > Q be a parameterization of the process 
(c, g )- The parameterization induces the mapping c|7: [0, a] -»c:si-> l{s)- The composition h o c\j 
is a function on [0, a]. The function h is increasing if and only if h o c|7 is increasing. Let h be in the 
ideal lo(c, go) C A(c, qo). For k < c, the jet ] k h G l fc o( c , Qo) is sa id to be positive if the jet j k (hoc\-f)(0) 
is positive. The jet } k h is said to be negative if the jet j k (h o c|7)(0) is negative. Being positive or 
negative is a property of the jet \ k h independent of the choice of the parameterization. The following 
proposition is an adaptation of Proposition 2. 

Proposition 3. Let k < c. Ifajet} k h in the ideal l fe o(c, q ) is positive, then the germ \ c h is increasing. 
If the jet \ k h is negative, then the germ \ c h is decreasing. 

The proposition implies that if the germ \ l h of a function h is increasing, then the jet \ k h is non 
negative for each k < c. 

Let g be a configuration in the set C° an let k G N U {00}. If for each jet j k (c, q ) of an admissible 
process (c, go) the jet j fe W( c , go ) is positive, then the germ j c W( c _ qo y Hence, go is a stable equilibrium 
configuration. We have obtained a sufficient condition for a point g G C° to be a stable local 
equilibrium configuration for each k G N U {00}. 



7 



If g is a stable equilibrium configuration, then the germ j c W( c ,<j ) 1S increasing for each germ j c (c, g ) 
of an admissible process (c, go)- It follows that for each k G N U {00} and each jet j k {c 1 go) of 
an admissible process (c, go) the jet j fe W( c ,g ) 1S non negative. This results in a series of necessary 
conditions for a point g G C° to be a stable local equilibrium configuration. Neither of these conditions 
is sufficient. Variational principles of classical physics are based on the necessary equilibrium condition 
of order k = 1. 

6. Composition of static systems and control. 

Static systems can be composed if they share the same control configuration space. The equality 
of configuration spaces is usually the result of a suitable choice. Let two systems be characterized by 
work functions W\: Ci — > R and W%: C2 — > R defined on sets Ci and C2 of admissible processes in the 
same configuration manifold Q. The set C = Ci fl C2 is the set of admissible processes of the composed 
system and the work is the function W — Wi\C + \V2\C. Coupling the system to other systems with 
the same configuration space is a form of control. The function W is used to find the equilibrium 
configurations of the controlled system. 

Composition of constrained systems may encounter difficulties. Examples will be given. 

7. Partial control. 

We have considered control of static systems through interaction with systems with the same con- 
figuration space. This is not always the case. One can in general associate three distinct configuration 
spaces with a static system: the internal configuration space Q, the control configuration space Q, and 
the observed configuration space Q. There are differential relations connecting the three spaces. We 
will consider the cases when a static system with a configuration space Q is controlled by external 
devices in a configuration space Q. The relation between the two spaces is a differential fibration 
n:Q^Q. The configuration space Q of the controlled system is the internal configuration space and 
the configuration space Q of the controlling devices is the control configuration space. We will refer 
to such situations as cases of partial control. The observed configuration space Q will coincide either 
with Q or with Q. 

Let a controlled system be represented by a work function W: C — > R defined on a set of admissible 
processes C c PQ and let a controlling device be represented by a work function W: C — > R defined 
on a set of admissible processes C C PQ. An admissible process of the combined system is a process 
c G C such that c' = 77(c) G C. Let C G C be the set of such processes. The work of the combined 
system is the function 

F:C^l:c^ W(c) + W'{r){c)). (34) 

We denote by PvQ the set of processes (c, q ) such that 77(c) is a single point set. Let CV be the 
intersection C (~l PvQ- We introduce the concept of the critical set. A configuration g G C° is in the 
critical set Cr(W, 77) if for each germ j c (c, g ) G Cy\ Q the germ j c u>( Ci<?0 ) is increasing. 

8. Realistic static systems. 

8.1. Constraints determined by differential equations. 

From the set C G PQ of admissible processes of a static system we extract the sequence of sets 

C° G Q, C 1 G P 1 ^, . . . , C' k G P fc Q, ... , C°° c P°°Q, C c G P C Q. (35) 

The fc-jet of a process is in C k if it has a representative in C. The set C° was already introduced 
earlier. Relations 

n k ' kQ (C k ) = C k ' (36) 

hold for k' ^ k and 

7 r kQ (C k )=C° (37) 

for k > 0. 
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Let (c, g ) be a process with the property that for each configuration q G (c, q ) with the exclusion 
of the terminal configuration qi the jet j fe (c| g , q) of the process (c| g , q) based on the arc 

C| g = {g'Gc; q'^q}. (38) 

is in C k . We denote by C k the class of all such processes. If a process (c, go) with terminal point gi is 
in C fe , then j fe (c| g , g) G C k for each g G c \ {gi}. It follows that j fe '(c| g ,g) = ir k ' k Q(j fe (c| g , g)) G C k ' for 
each g G c \ {gi}. Hence, (c, g ) G C k . We have established the inclusion 

C k ' D C k (39) 

for k! sC k. 

If (c, g ) is an admissible process, then for each configuration q G (c, go) with the exclusion of the 
terminal configuration gi the process (ci g ,g) is admissible and the germ j c (c| g , g) is in C c . It follows 
that C G C c . Let the germ j c (c| g ,g) be in C c for each configuration q G c \ {gi}. We show that the 
process (c q « | , g ) is admissible for each configuration g" in the interior of c. For each g G c \ {gi} there 
is a configuration q' G (c, g ) such that g' > g and the process (ci gg /i,g) based on the arc 

C| gg ,| = {g" £c; ? ') q" ^ g} (40) 

is admissible. This is a consequence of j c (c| g ,g) G C c . The union 

U ( c to'l><?) (41) 
<?ec\{ 9 i} 

is an admissible process and for each g" in the interior of c the process (c g j|,go) is a subprocess 
of the union. Hence, (c g «|,go) is admissible. The process (c, g ) itself is admissible since it can be 
approximated by admissible processes. We have proved the inclusion C c C C. 

It may happen that the equality C = C k holds for k < c. If it holds for k and not for k — 1, then the 
set C is said to represent constraints of order k. The set C k is interpreted as a differential equation of 
order k. Admissible processes are obtained by solving this equation. Constraints of order are said 
to be holonomic. Constraints of order 1 are said to be non holonomic. Constraints of higher order are 
not usually discussed since their presence is not apparent if only first order criteria of equilibrium are 
considered. 

8.2. Jets of processes with volume, the work integral. 

A k-jet of a process with volume is a pair (j fe (c, g ), v), where v is a vector belonging to the orientation 
o of (c, go) at go- We denote by V fe Q the space of k-jets with volume. For k < c it make sense to 
associate A;-jets with volume with the terminal configuration of a process. Spaces V°Q and V 1 are 

o 

identified with the bundle TQ C TQ of tangent vectors with the zero vectors removed. This is possible 
since in the pair (j°(c, go),f) or (j 1 (c, go),f) the vector v already contains the complete information 
about the jet. We denote by V k the fc-jets with volume associated with admissible processes. Note 
that if (j fe (c, g ), v) G V k , then (j fe (c, g ), \v) G V k for each A > 0. Sets V° and V 1 are interpreted as 

o 

a subset V C TQ. 

The work of a realistic static system is defined in terms of a work form 

tf: V k -» M. (42) 

The work form is positive homogeneous in its vector argument: if A > 0, then $(j fc (c, g ), \v) = 
Xi3(i k (c, g ), v). The work form is used to define the work function W(c, 9o ) along an admissible processes. 
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Let gi be the terminal point of the process (c, q ). For each configuration q £ c \ {gi} the value of the 
work function is the integral 

">(«,«,)(</) = I <? (43) 

J (c\q,q ) 

of the work form defined as the Ricmann integral 

o 



f 

Jo 



tf(j fc (c| 7(s) , 7 ( S )),t 7 ( S ))d S (44) 



in terms of a parameterization 7: [0, a] — ► Q. Homogeneity of 1? makes the integral independent of the 
parameterization. The work of the process is the limit 

W(c,q ) = lim w {c , qa) (q). (45) 

The work form of a system of the usually considered type is a positive homogeneous function 

1?: V -> R. (46) 

This corresponds to k = or k = 1. 

9. The principle of virtual work. 

Variational principles of classical physics are versions of what is known in statics as the principle of 
virtual work. This principle consists in applying the first necessary equilibrium condition to a static 
system. It is usually assumed that the constraints are of order or 1 and that the work is represented 
by a work form 

0: V -> R. (47) 

Elements of V are the admissible virtual displacements. The principle of virtual work states that the 
inequality 

d(v) > (48) 

holds for each virtual displacement v tangent to an admissible process initiating at a configuration of 
equilibrium. 

The principle of virtual work can be extended to more general realistic systems. If the system has 
constraints of order I and the work is defined in terms of a work form 

1?: V k R, (49) 

then the principle of virtual work (the first order necessary condition of equilibrium) states that the 
inequality 

#(j k (c,q o ),v)>0 (50) 

holds for each fc-jet with volume (j fc (c, q ), v) belonging to an admissible process (c, go) initiating at a 
configuration of equilibrium. 

10. Potential systems and the Legendre transformation. 

A potential order 1 of a potential static system is a diffcrcntiablc function U on the configuration 
space Q. A potential system is unconstrained. For each process (c, q ) e C = PQ with terminal point 
qi the function W( c , 90 ) i s defined by 

w {cm) {q) = U{q)-U{q Q ) (51) 
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for q 6 c \ {qi} and the work W(c, qo) = W( Ciqo )(qi) is the limit 



W(c,q ) = lim w (c . 9o) ((?). 



(52) 



The set F of a potential system is the tangent bundle TQ and the function -d is the function 

z9: TQ — > 18: j) i — ► (du, u) (53) 

constructed from the differential df/ of the potential. The principle of virtual work for a potential 
system states that the equality 

{dU,v)=0 (54) 

holds for every virtual displacement from a configuration of equilibrium. The appearance of the equality 
in place of the inequality is due to the reversibility of virtual displacements: if v is a admissible virtual 
displacement from a configuration of equilibrium, then so is — v and (dll, —v) = — {dU, v) 

o 

The Legendre transformation associates with a static system represented by a set V C TQ and a 
function fl: V — ► R the constitutive set S of the system defined as 

S ={f e T*Q; ttq(/) e C°,t?(t>) - (./» ^ for each u e V such that tq(v) = ttq(/)} (55) 

The set C° is obtained by applying the tangent projection tq to virtual displacements. If constraints 

o o 

are holonomic, then V = TC°. A virtual displacement is in TC° if there is a curve £:M — ► Q and a 
number e > such that £([0,e]) C C° and w = t£(0) 7^ 0. The inequality in the definition of S is 
replaced by an equality in the case of reversibility. 

The physics of the Legendre transformation is that of control of a static system by means of potential 
external devices. The covector / = — dU(wQ(f)) is the external force applied to the controlled system 
by the potential device with potential U. The constitutive set S characterizes the response of a static 
system to control by potential devices only. Under certain conditions (convexity) this characterization 
is complete since the set V and the function d can by reconstructed from the constitutive set by the 
inverse Legendre transformation. 



11. Examples of static systems and Legendre transformations. 

Configuration spaces of most systems considered here are affine spaces. If Q is an affine space 
modeled on a vector space W, then the tangent bundle TQ is identified with the product Q x W, the 
cotangent bundle is identified with Q x W* and the canonical pairing is the mapping 

( , ): (Q x W*) x (Q x W): ((q, /), (q, w)) ~ (/, «;>. (56) 

We denote by q\ — qo the vector associated with the points qo and q\ . 

For the sake of simplicity we will not exclude zero vectors in the set V of virtual displacements. 

Example 1. A material point with configuration q in the affine Euclidean space Q of Newtonian 
mechanics is tied to a fixed point qo with a spring of spring constant k. This is an unconstrained 
potential system. The internal energy is the function 

U:Q^R:q^^\\q-qo\\ 2 . (57) 

The set 

S={(qJ)£QxW* ] f = kg(q- q )} (58) 

is the constitutive set. ▲ 
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Example 2. Let go be a point in Q and let lo: W x W — > R be a bilinear mapping. The function 

x W -> R: (<?, w) uj(q - q ,w) (59) 
used in the principle of virtual work produces the constitutive set 

S = {(q,.f)eT*Q; w(q-q ,v)-(f,v)=0 for each v e W} 

= {(q,f) eT*Q; w(g - g , •) - / = 0} . (60) 
The function $ is obtained from the differential of the quadratic function 

U:Q^R:q\-^^u>(q-q ,q-qo) (61) 

only if u is symmetric. ▲ 

Example 3. Let Q be the configuration manifold of a static mechanical system and let p: TQ — > T*Q 
be a mapping with the properties of a Euclidean metric. This means that p defines an isomorphism of 
vector bundles 



(62) 



the symmetry relation 

(p(u),v) = (p(v),u) (63) 

holds for any pair (u, v) e TQ x TQ, the inequality 

( t q, t q) 

(p(v),v)>0 (64) 

holds for each »^0. The function 

§:TQ^R:v^ y/(p(v),v) (65) 

is positive homogeneous on fibres of the tangent fibration. It represents the virtual work of a virtual 
displacement v due to friction. The principle of virtual work is the inequality 

y/{p(v),v)-(f,v)>0. (66) 
It is satisfied if / is in the constitutive set S and (/, v) eT*Q x TQ. 

Let / 6 T*Q be in the constitutive set. By using v = p^ 1 (,f ) in the principle of virtual work we 
arrive at the inequality 

{f,p-\f))^y/(f,p-Hf))- (67) 

Hence, 

(68) 

The inequality 

(f,v) < y/{f,p- l tf))y/(p(v),v). (69) 
is the result of the Schwarz inequality 



( P (U),V) < y/(p(u),u)y/{p(v),v) (70) 

applied to the pair of vectors u = p~ 1 (f) and v. If (/, p~ 1 (f)) ^ 1, then 

(f,v) < y/(p(v),v). (71) 

Hence, / € S. 

The constitutive set of the system is the set 

S = {feT*Q;{f, P -Hf))<l}- (72) 

▲ 
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Example 4. Let a material point with configuration q in the Euclidean afhne space Q be tied with a 
rigid rod of length a to a point with configuration qo . The configuration q is constrained to the sphere 

C° = {qeQ; \\q-q \\=a}. (73) 

This is a system with holonomic bilateral constraints. The set 

V = {(q,5q)eQxW; ||«-«b||=o, (g(q - q ), Sq) = 0} (74) 

of admissible virtual displacements is the tangent set TC° of the holonomic constraint C . With the 
work form ■& — the constitutive set is the set 

S = {(?,/) EQxW*; \\q-q \\=a, f = a- 2 (f,q-q )g(q-qo)}. (75) 

▲ 

Example 5. Let g be a point in the Euclidean affine space Q of dimension 3 and let u\, u 2 and u 3 
be mutually orthogonal unit vectors in the model space W. Let the configuration q e Q of a material 
point be constrained to the set 

C° = {q G Q; (ff(ui), g - g > ^ 0, 0/(« 2 ), g - g ) > 0} . (76) 
This is a system with holonomic one-sided constraints with the tangent set 
V = {(q, Sq)eQx W; (g( Ul ),q - q ) > 0, (g{u 2 ), q - q ) > 0, 

0/(ui), <5g) > if 0/(«i), g - qo) = 0, (5(«2), Sq) > if ( 5 (w 2 ), q - q ) = 0} . (77) 
of C° as the set of admissible virtual displacements. With = the constitutive set is the set 
S= {(q,f) G Q x VK*; ^(ui), g - g > > 0, («/(« 2 ), g - q ) > 0, 

/ = if (g{ui),q - q ) > and (g{u 2 ),q - qo) > 0, 
(/,m)<0, (/,U2>=0, and (/,« 3 ) = 

if (g(ui),q -qo) =0 and (g{u 2 ),q - q ) > 0, 
(/,ui)=0, (/.ua) <0, and (/,« 3 ) = 

if {g{ui),q - qo) > and (g(u 2 ),q - q Q ) = 0, 
(j>i) <0, (/,u 2 > <0, and (/,« 3 ) = 

if (g{ui),q-qo) =0 and (g{u 2 ),q-q ) = 0} . (78) 

▲ 

Example 6. Let Q = X x U be the configuration space of a skate. The space X is an affine plane 
modeled on a vector space W and D is the projective space of directions in X. We use a Euclidean 
metric in X represented by a mapping g:W — > to identify the space D with the unit circle 

D = { V eW; (g(<p), V ) = l}. (79) 

This is the only use we make of the metric. The tangent bundle TQ is identified with X x W x TD, 
where 

TD = {(if, 5<p)eDx W; {g((p),6<p) = 0} . (80) 

The skate is a system with non holonomic constraints. The set C° is the entire space Q. The constraint 
consists in restricting virtual displacements in X to those parallel to the direction specified by an 
element of D. Thus 

V = { (x, Sx, if, Sip) e X xW x TD; 3 km Sx = kip} . (81) 

The cotangent bundle T*Q is the product X x W* x T*D. For each ip G D the fibre T V D is a subset 
of W. Hence, the cotangent bundle T*D can be specified as the set of pairs (ip,r), where ip G D and 
t is in the quotient space W* /T°D. The set 

S = { (x, f, ip, t) G T*Q; (/, + (t, Sip) = for each (x, Sx, ip, Sip) G V} 

= { (x, /, ip, t) G T*Q; (/, ip) = 0, r = 0} (82) 
is the constitutive set of the system. ▲ 
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Example 7. Let Q be the Euclidean affine space of Newtonian mechanics. The model space for 
Q is a vector space W of dimension 3. The Euclidean structure is represented by a metric tensor 
g:W^W*. 

The example gives a formal description of experiments performed by Coulomb in his study of static 
friction. Let a material point be constrained to the set 

C° = {qeQ: (g(k),q-q ) ^0}, (83) 

where qo is a point in Q and k £ W is a unit vector. The boundary 

dC° = {qeQ; (g(k),q-q )=0} (84) 

is a plane passing through qo and orthogonal to k. In its displacements on the boundary the point 
encounters friction proportional to the component of the external force pressing the point against 
the boundary. The system is characterized by the virtual work function ■§ = defined on the non 
holonomic constraints 

V = [(q,5q) GTQ; (g(k), q - q ) > 0, (g(k),6q)) > vy/\\5q\\* - (g(k),6q))* 

if (g(k),q-q ) = o}. 

(85) 

The principle of virtual work states that (q, f) is in the constitutive set S if and only if the inequality 

</,<*?) <0 (86) 

is satisfied for each (q, Sq) £ V. 

If (g(k), q — qo) > 0, then a pair (q, /) £ T*Q is in the constitutive set S if and only if / = 0. 

We consider pairs (q, /) with (g(k),q — qo) = 0. If / = — ||/||fc, then (q, /) is in the constitutive set 
and \\f\\ 2 - (/, k) 2 = 0. Let (q, f) be in the constitutive set and let ||/|| 2 - (f,k) 2 ^ 0. The virtual 
displacement (q,Sq) with 

Sq = g~\f) ~ </, k)k + zVll/ll 2 - (/, k) 2 k (87) 

is in V since 

(g(k),5q) = ^\\f\\ 2 -(f,k) 2 . (88) 

From the principle of virtual work and 

(/, Sq) = H/ll 2 - (/, k) 2 + y/\\f\\ 2 -(f,k) 2 (f, k) (89) 
it follows that 

ll/ll 2 - (/, k) 2 + VU¥~ifl¥(f, k) < o (90) 

and 

y/\\f\\ 2 -(f,m + "(f,k)^0 (91) 

since H/ll 2 - (f,k) 2 ) >0. 
The Schwarz inequality 

- (. 9 (fc), U )(. 9 (fc)^} < ^/\\u\\ 2 - (g(k),u) 2 ^\\v\\ 2 - (g(k),v) 2 (92) 

for the bilinear symmetric form 

(u,v) i ^ {g(u),v) - (g(k),u)(g(k),v) (93) 
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applied to the pair (g 1 (f),Sq) leads to the inequality 

(/, Sq) - (/, Sq)(g(k),Sq) < y/\\f\\* - (f,k)*y/\\6q\\2 - (g(k),6q) 2 - (94) 

If (g(k),q- 90) = 0, Vll/P-a m + v(.f, k) < 0, and (g(k),6q)) > vy/\\6q\\* - {g(k),6q))*, then 
(/> Sq) ^ 0. Hence, (q, /) is in the constitutive set S. 
We have shown that the set 

S={(q,f)eT*Q; (g(k), q - q ) > 0, f = if (<?(£;), g - q ) > 

and V\\f\\ 2 -(f, k) 2 ) + u{f, k) < if (g(fc), g - 90} = o} (95) 

is the constitutive set of the system. 

▲ 

12. The Legendre transformation for partially controlled systems. 

Let Q be the internal configuration space of a system, let Q be the control space and let rj:Q~^Q 
be a differential fibration relating the two spaces. The work function V — > K is defined on a set 
V C TQ. The system is partially controlled by potential systems. If Q = Q, then 

S =((«,/) €Q x T*Q; g e Tq(V), - (f,Tr](v)) > 

for each v such that Tq(u) = g| . (96) 

If Q = Q, then the constitutive set S C T*Q is the image pr T *Q(S) of S by the canonical projection 

pr T * Q :Q x T*Q -> T*Q. (97) 

The concept of the critical set adapted to first order equilibrium criteria is useful. The critical set 
is the set 

Cr = {9 e Q; g G 7-q(V), > for each v eV such that Tq(w) = q and = o| . (98) 

13. Examples of the Legendre transformation for partially controlled systems. 

We give examples of partially controlled systems in the affine space Q of Newtonian mechanics 
modeled on a Euclidean vector space W with a metric tensor g:W — > W* . The tangent bundle TQ is 
identified with the product Q x W and the cotangent bundle is identified with Q x W* . The tangent 
and the cotangent bundle of the product Q x Q are identified with QxQxWxW and QxQx W* x W* 
respectively. 

Example 8. Three material points with configurations q , qi, and q 2 in the affine space Q are 
interconnected with springs with spring constants k Wl fc 2 o, and fc 2 i- The virtual work of the system is 
derived from the internal energy 

U: Q x Q x Q - R: (90,91,92) ~ ^fhi - 9o|| 2 + ^||<?2 - 9o|| 2 + ^||?2 - 9i|| 2 (99) 

of the springs. The point 90 is fixed and not controlled. The configuration space is the product Q x Q. 
The two points 91 and 92 are not constrained. The constitutive set is the set 

S = {(9i,<72,/i,/2) e T*(Q x Q); fi = k 10 g{qi -90) - k 21 g(q 2 -91), 

h = feo.9(<?2 - 9o) + ^215(92 - 9i)} • (100) 
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If the configuration of q 2 is not controlled but observed, then Q x Q is the internal configuration 
space, Q is the control configuration space and the canonical projection 

wv-Q x Q -> Q- (91,92) >-»• 9i (101) 

is the relation between the two spaces. The critical set 

Cr = {(q u q 2 ) E Q X Q; k 20 g(q 2 - q ) + k 21 g(q 2 - ?1 ) = 0} (102) 

is the image of the section 

fi:Q —> Q X Q:q>-> (q,q + k 2 i(k 20 + fc2i)~ 1 (g - q )) (103) 

of the projection pr\. The set 

S = ((9i,92,/i) € Q x Q x W*; q 2 = q + ^ (gi - <? ), 
I fc 20 + K21 



A -(*» + s£sr) ««-»»} (104 » 



is the constitutive set. 

If the configuration of q 2 is neither controlled nor observed, then the set 

S = G Q x /1 = + ^J^ ) .9(91 - 9o)| (105) 

is the constitutive set. Note that the presence of the material point with configuration q 2 can be 
ignored. This is due to the fact that the critical set is the image of a section of the projection (101). A 

Example 9. The present example gives a simplified discrete model of the buckling of a rod. One 
end of the rod is a point in the affine space Q with configuration qi constrained to the line 

L = {9 € Q; q - qo = (g{u),q - qo)u} (106) 

through a point go in the direction of a unit vector u. The other end is a point with configuration q 2 
constrained to the plane 

P={qeQ; (g(u),q-q )=0} (107) 

through qo perpendicular to u. The rod can be compressed or extended in length but not bent. Its 
relaxed length is a and the elastic constant is k. The buckling of the rod is simulated by displacements 
of its end point in the plane P tied elastically to the point q with a spring of spring constant k! . The 
configuration space is the product Q x Q with holonomic constraints represented by 

C° = {(91,92) EQxQ; q 1 -q Q = (g(u) , q 1 - q )u, (g{u),q 2 -q ) =0}. (108) 

The set 

V = {{qi,q 2 ,8q 1 ,5q 2 ) £ T(Q X Q); q - q = {g{u),qi - q )u, 

(g(u),q 2 - 9o) = 0, 5qi = {g{u),Sqi)u, {g{u),Sq 2 ) = } (109) 

of admissible virtual displacements is the tangent set of C°. The internal energy of the system is the 
function 

U:C° - R: {q u q 2 ) ~ - q 2 \\ - af + ^\\q 2 - q Q \\. (110) 
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This expression for the internal energy excludes certain non realistic configurations. The set 
S = {(<7i,<?2,/i,./2) G T*(Q x Q); qi-q = (g(u),qi - q )u, (g(u),q 2 - q ) = 0, 

(fi,u) = k (l - a||g 2 - gilT 1 ) (g{qi - go),"}, 

f 2 - (f 2 ,u)g(u) = (k + k'- ka\\q 2 - q 1 \\- 1 ) g(q 2 - q )} (111) 

is the constitutive set of the system. 

If the configuration q 2 is not controlled, then the critical set is the union of sets 

Cry = {(qi,q 2 ) G Q x Q; qi - q = (g(u),qi - q )u, q 2 = q } (112) 

and 

Cr 2 = {(qi,q 2 ) G Q x Q; q 1 - q = {g(u),qi - q )u, 

(g(u),q 2 -q )=0, (k + k')\\q 2 - qi \\ = ka}. (113) 

If the configuration q 2 is not controlled but observed, then the constitutive set is the union of sets 

Si = {(qi,q2,.fi) eQ xQ x W*; q 1 - q = (ff(u),gi - q )u, 

q2=q\, (fi,u) = k (l - a||gi - golT 1 ) {g(qi - qo),u)} (H4) 

and 

S 2 = {(gi,(?2,/i) e Q x Q x W*; qi-q = (g{u),qi - q )u, (g(u),q 2 - g ) = 0, 

(k + k')\\q 2 - qi \\=ka, (fuu) = -k'(g( qi -q ),u)}. (115) 

If the configuration q 2 is neither controlled nor observed, then the constitutive set is the union of 
sets 

Si = {(<?i,/i) G Q x W*: q 1 -q = (g(u),qi - q )u, 

(f u u) =k(l- a\\ qi - qoir 1 ) (g(qi - q ),u)} (116) 

and 

= G Q x W*; qi-q = (g(u),qi - q )u, 

(k + k')\\ gi - q \\ < ha, (fi,u) = -k'(g( qi -q ),u)}. (117) 

The constitutive set S — Si U S 2 is not a submanifold of T*Q. ▲ 

Example 10. A material point with configuration q 2 in the affine space Q is connected to a fixed 
point g with a rigid rod of length a. A second material point with configuration q\ is tied elastically 
to q 2 with a spring of spring constant k. The configuration space is the product QxQ with holonomic 
constraints represented by 

C° = {(qi,q 2 )eQxQ; ||<ft - (toll = a}. (118) 

The set 

K = {(«i, (ft, *gi,<fcfc)eT(QxQ); \\q 2 - q \\ = a, (g(q 2 - Sq 2 ) = 0} (119) 
is the tangent set of C°. The set 

S = {(<?i, 92, /i,/2) G T*(Q x Q); \\q 2 -g 1| = a, /i = fc.g(gi - g 2 ), 

/2 - fc5(?2 - <?i) = a~ 2 (/2 - fcgfe - gi), qi - qa)g(q2 - qa)} ■ 

(120) 
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is the constitutive set and 

Cr = {(q 1 ,q 2 )eQxQ; \\q 2 - q Q \\ = a,\\ qi - q a \\(q 2 - q a ) = ±a( qi - q )} . (121) 

is the critical set. 

If the configuration q 2 is not controlled but observed, then 

S={( qi ,q 2 ,fi) eQxQxW*; \\q 2 - (toll = a, 

hi - 9o!l(<?2 - go) = ±a(gi - go), A = Mgi - 12)} ■ (122) 

is the constitutive set. 

If q 2 is not observed, then we have the constitutive set 

S={(gi,/i)eT*Q; UMI = ka if q x = q , 

h = k(l± a\\ qi - goir 1 ) g( qi - g ) if gi + g } . (123) 

This set is the image of the injective mapping 

p: K x R -> Q x W*: (w, r) ^ (q + rw, k(r - a)g(w)), (124) 
where K is the unit sphere 

K = {w e W; |HI = !}• (125) 
The space T(K x R) is identified with the set 

{(w, r, Sw, Sr) e W x R x (W © R); ||w|| = 1, (<?(w), <5w) = 0} (126) 

and the space T(Q x W*) is identified with Q x W* x Q x W* . The tangent mapping 



Tip: T(K xl)-t T(Q x VF*) 

: (w, r, <5u>, <5r) i— ► (<jo + ™, £;(r — a)g(w), Srw + rSw, kSrg(w) + k(r — a)g(Sw)) 



(127) 



maps injectively each tangent space of K x R at (w, r) into the tangent space of Q x W* at <p(w, r). 
It follows that ip is an injective immersion. 

The space TT*Q x TT*Q is identified with Q x W* x (W ®W*) x (W ®W*). The symplectic 

T*Q 

form on Q x W* is introduced as the bilinear mapping 

uq-- Q x W* x (W © W*) x(W® W*) -» R 

: (q, f, 5 iq , <5 2 g, <5 2 /) -> <^/, <5 29 ) - <<5 2 /, fcg). (128) 

The product T(K x R) x T(lf x R) is the set 

K xR 

{(w, r, <5iw, <5ir, 5 2 w, 5 2 r) e W x R x (VK © R) x (W © R); \\w\\ = 1, 

(g(w), S lW ) = 0, (g(w), S 2 w) = 0} (129) 

and the pull back form 

tp*w Q :T(K xR) x J(K xR)^l 

KxR 

: (w, r, Siw, Sir, S 2 w, S 2 r) i— > (k5\rg{w) + k(r — a)g(<5iw), <5 2 rui + r5 2 w) 

— (kS 2 rg(w) + k(r — a)g(8 2 w) 1 8irw + rSiiv) (130) 
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is zero. We have shown that the image S of (p is an immersed Lagrangian submanifold of T*Q. It can 
be shown to be embedded. The mapping 



7Tq o ip: K x M — > Q: r) g + rw (131) 

is not injective. The tangent mapping 

T(ttq o <p): TfTf xR)^ TQ: (w, r, Sw, Sr) (g + rw, drw + rSw) (132) 

is of rank 3 if r ^ and of rank 1 if r = 0. This indicates the presence of a Lagrangian singularity of 
S above the point go- 

Note that the critical set is not the image of a section of n. For each control configuration (x, y) we 
have two different internal equilibrium configurations if x 2 + y 2 > and an infinity of internal equilib- 
rium configurations if x 2 + y 2 = 0. The external force necessary to maintain the control configuration 
(x, y) depends on the internal configuration. Thus even if the internal configuration is not directly 
observed its presence can not be completely ignored. A 

14. The Legendre transformation of composed static systems. 

Let two static systems in a configuration manifold Q be characterized by sets (C°i, V\, Si) and 
(C°2, V2, 1S2S2) respectively. Let the composition of the systems be characterized by (C°, V, S). The 
constraints C° and V are the intersections C° = C°i n C° 2 and V = Vi D V 2 . The form fr. V -> R is 
the sum 1? = ^\V + d 2 \V. Sets S lq = Si n T*Q, S 2(J = S 2 n T*Q, and 5, = Sn T*Q are defined by 

Si 9 = {/eT g Q; = t?ij, (133) 

S 29 = {/ e T g Q; f\V 2g = t? 2 J , (134) 

and 

S,={/£T,Q; f\V q = d q }. (135) 

If /1 6 Si, and / 2 G S 2g , then /1 + / 2 G S g . This establishes the inclusion S\ q + S 2g C S 9 . Assuming 
that sets Vi g = V\ fl T g Q, = fl T g Q, and 1^ = 1^1 T 9 Q are subspaces of T 9 Q we will show that 
the inclusion S q G Si q + S 2(J holds as well. We choose complements V[ q of V q in Vi q , V 2 ' q of V" g in V2 q , 
and Vb g of V\ q + V2 q in T 9 Q. Relations 

T q Q = Fi, + V 2 \ + V 09 = + V[ q + Vo g (136) 

^n^HO}, v 2 ',n^ = {o}, K ,n^i, = {o}, (137) 

and 

v 2q nv{ q = {0}, v{ q nv Qq = {0}, v Oq nv 2q = {0} (138) 

hold. If / G Sq, then / is the sum of covectors /1 G Si g Q an d f 2 £ S 29 Q characterized by 

A\Vi q = # lq , fi\Vl q = f\Vl q -0 2q \Vl q , h\V 0q = l -f\V Qq , (139) 

and 

f2\V 2q = #2 q , f2\V{ q = f\V 2 ' q -$ lq \V{ q , f 2 \V 0q = \f\V 0q . (140) 

We have established rules of composition of static systems is sufficiently simple cases. These rules 
are expressed by 

c° = c°i n c° 2 , (141) 
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v = v x nv 2 , 



(142) 



= 0i\V + & 2 \V, (143) 

and 

S q = Si, + S 2q (144) 

for each q G C°. 

Two submanifolds C\ and C 2 of a manifold Q are said to have clean intersection if the intersection 
C\ n C2 is a submanifold and TC\ (~l TC2 = T(Ci PI C2). Examples show that the simple rules listed 
above do not apply to the composition of systems with holonomoic constraints if the intersection of 
constraints is not clean. 

Example 11. We consider the composition of two holonomic systems of the type described in 
Example 4. The constraints and the constitutive set for the first system are represented by the sets 

C\ = {qeQ; \\q- qi \\=a}, (145) 

Vi = {(q,6q)€QxW; \\q - qi\\ = a, (g(q - qi), Sq) = 0} , (146) 

and 

51 = {(<?,/) eQxVK*; h-q2\\=a, f = a - 2 (f, q - qi )g(q- gi )} . (147) 

For the second system we have 

C° 2 = {qeQ; \\q-q 2 \\=a}, (148) 
V 2 = {(q,Sq)eQxW; \\q-q 2 \\=a, (g(q - q 2 ), Sq) = 0} , (149) 

and 

5 2 = {(?,/) eQxW*; \\q-q 2 \\=a, f = a~ 2 </, q - q 2 )g(q - q 2 )} . (150) 

If the distance \\q 2 — q\\\ between the centres of the spheres C°i and C° 2 is less than 2a, then the 
composed system is a system with holonomic constraints. The intersection 

C° = C\ n C\ = {qeQ; \\q - gi|| = a, \\q - q 2 \\ = a} (151) 

is clean. The intersection 

V = V 1 nV 2 = {(q,6q)eQxW; ||g-«i||=0, \\q - q 2 \\ = a, 

(g(q - qi ), Sq) = 0, (g(q - q 2 ),Sq) = 0} (152) 

is the tangent set TC°. The constitutive set 

S = {(«,/) eQxW*; \\q-qi\\=a, ||g-<6||=o, </,^}=0 for each 

Sq G W such that (g(q — q\), Sq) = and (g(q — q 2 ), Sq) = 0} . (153) 

is obtained from the principle of virtual work. At each q e C° the set 

S q = {feW*; (q,f) 6S} (154) 

is the sum 

{/ e W* ; ( q j) e Si} + {/ e W* ; (<?,/) e S 2 } . (155) 

The simple rules are followed in this case. 
If \\q 2 — qi || = 2a, then the set 

C° = C 1 nC° 2 = {qeQ; ||«-gi||=a, ||g-ga|| = a} (156) 
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has only one element q = q\ + \{q 2 — qi). The intersection V\ n V 2 is the set 

l {q, Sq)eQxW; q = qi + ^(q 2 - gi), (g(q 2 - qi ),5q) = j . (157) 

This is not the tangent set of C°. The intersection of C°i with C° 2 is not clean. With 

V = TC° = {(q,f)&QxW*; q = q, + \{q 2 - q x ), Sq = o| (158) 

the principle of virtual work produces the constitutive set 

S = {(<?,/) G Q x FF*; g = gi + ife - ft) J • (159) 

This is not the correct constitutive set for the composed system. The reason of this failure is that the 
approximative assumption of perfect rigidity of the separate constraints is no longer realistic in the 
case of composition which is not clean. To obtain a complete description of the composed system the 
precise elastic properties of the constraints must be known. A partial characterization of the system 
is provided by the constitutive set 



S = /) € Q x W*; q = Ql + ^(q 2 - qi ), f = a- 2 (f, q - qi )g(q - 



(160) 



generated by the non holonomic constraints V = V\ fl V 2 ■ Note that these constraints are not integrable 
since the inclusion V C TC° does not hold. The simple composition rules can not be followed. ▲ 

15. Configuration spaces for the dynamics of autonomous machanical systems. 

We consider the space Q(M|M; [to, ti]) of differentiable mappings 7: [to, t\] — > R of an interval [to, t{\ c 
R in a manifold M. This space serves as the configuration space for dynamics in the time interval 
[tojti]- It is not a differential manifold. Basic constructions of differential geometry are nevertheless 
possible. Tangent vectors and covectors can be defined if a class of functions and a class of curves are 
chosen to act as differentiable mappings. 

From a function 

ip: TM — > R (161) 

we derive the function 

rti 

F:Q(M|R: [t ,h]) -> R:7 ^ / ip o t 7 . (162) 

Jto 

Functions on Q(M|R; [to,ii]) obtained by this construction are declared differentiable. 
Starting with any differentiable mapping 

R x [t , t\] — > M (163) 

we introduce the curve 

r: R - Q(M|R; [t , h]): s h-» x («, ■ ) (164) 

in Q(M|R; [to,*i]) considered differentiable. 

The restriction of a differentiable curve to an interval [0, a] is a parametrization of a process. Virtual 
displacements are tangent vectors defined as equivalence classes of differentiable curves. Curves T and 
r' are equivalent if 

(F o r')(0) = (F o r)(0) (165) 
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and 

D(For')(0) =D(For)(0) (166) 

for each difTerentiable function F. Let T and V be derived from mappings x : ^ x [to,ti] — > M and 
x':R x [t , ^i] — > M and let F be constructed from a mapping (f:TM — > R. Equivalence conditions 
involve the function 

FoT:R -> R:s ^ / <p o tx(s, • ) (167) 
J to 

and a similar function with x'- Equivalence conditions will be satisfied if x an d x' produce the same 
mapping 

w:[to,h] -f TM:ti-f tx(-,t)(0). (168) 

This mapping is a natural choice of a representative of the equivalence class of T. 

A covector is an equivalence class of pairs (q,F) of a configuration q G Q(M|R; [to, fi]) and a 
difTerentiable function F on Q(M|R; [t ,ti]). Pairs (q,F) and (q',F') are equivalent if = g and 

D(F' o T)(0) = D(F o T)(0) (169) 

for each difTerentiable curve T. 

Let F, F' and L be derived from mappings tp:TM -> R, y?: TM -> R and x:R x [t ,ti] -> M 
respectively. It follows from the formula 



D(FoT)(0) = / (d(f,tw) (170) 
Jto 



that the mapping dip o tx(0, •) could be used to represent the equivalence class of (q,F). This is 
not a convenient representation since tw is not a generic mapping from [t , ti] to TTM. A better 
representation is derived from the equality 

D(For)(o) = - / '(^ot^fo, -)^) + ((^°tx(o, -)).^i)>-((^°tx(o, -))Mh)), (171) 

where 

£,tp:T 2 M -^T*M (172) 

and 

Ttp-.TM ^>T*M (173) 
are mappings defined in [DA] . The equality suggests that a covector should be represented by a mapping 

/:[t ,ti]->T*M (174) 

and two covectors po and p\ in T*M. The pairing between a covector represented by {f,Po,Pi) and a 
vector represented by w is expressed by 

((/,Po,Pi)> = - / \f,w) + <Pi,u>(*i)> - <Po,M*o)>. (175) 

Jt 

In dynamics the term action replaces the term work in statics. For simple autonomous systems 
usually considered in analytical mechanics all processes are admissible and the action of a process is 
the difference of values at the end points of the process of a function 



rh 

L:Q(M|R; [t ,ti]) -» R: 7 ^ / A o t 7 (176) 

Jtn 



/to 
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derived from a Lagrangian 

A: TM — > R. (177) 

The variational principle 

/ 1 (dA, tw) = - [ 1 (/, w) + (pi, w(ti)) - (po, u>(to)} (178) 
produces the Euler- Lagrange equations 

£Aot 2 7 = / (179) 
in [iojti] and the momentum- velocity relations 

(TAot 7 )(t )=Po (180) 

and 

(?Ao7)(ti)=pi. (181) 

The formulation of analytical mechanics in a finite time interval has its infinitesimal limit with the 
interval shrinking to a point. The passage to the infinitesimal limit is greatly simplified by observing 
that in a slightly modified version of the formulation the interval [t , t\] is used exclusively as a domain 
of integration and can be treated as a current or distribution. Replacing the interval by the Dirac delta 
St at t e R produces the infinitesimal limit. We will describe the modified version of the formulation 
using currents in R in place of intrvals. 

We introduce an equivalence relation in the set of pairs (c, 7) , where c is a current in R and 7 is a 
curve in M. Pairs (0,7) and (c',7') are equivalent if c' = c and 

ip o fy' = / ip o t 7 (182) 



for each differentiable function p: TM — > R. The space Q(M|R; c) of equivalence classes is the config- 
uration space corresponding to the current c. The equivalence class of a pair (c, 7) will be denoted by 

q(c,7)- 

Each differentiable function ip: TM — > R defines a function 

F:Q(M|R;c) R:q(c, 7 ) ^ J ipot-f (183) 

considered differentiable. 
From a mapping 

XtRxR^M (184) 

we construct the curve 

r: R -> Q(M|R; c): s ^ q(c, x(s, ■ )) (185) 

in Q(M|R;c) differentiable by declaration. 

Tangent vectors are defined as equivalence classes of differentiable curves. Equalities (165) and 
(166) define equivalence of curves. The space of tangent vectors is denoted by TQ(M|R;c). A natural 
representative of a tangent vector is an equivalence class of pairs (c, w), where c is a current in R and 
w is a curve in TM . Pairs (c, w) and (c', w') are equivalent if c' = c and 

/ (d<p,tw') = I (d<p,tw). (186) 

The space of equivalence classes will be denoted by Q(TM|R; c). 
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It is convenient to consider the construction of spaces Q(M|R; c) and Q(TM|R; c) as the results of 
applications of a covariant functor Q( • |R; c) to manifolds M and TM. This functor applied to the 
morphism tm ■ TM — > M results in the morphism 

Q(t m |K;c):Q(TM|R;c) Q(M|R; c): q(c, w) ^ q(c,r M o w) (187) 

which can be treated as a vector fibration. The morphism 

tq(m|R;c) : TQ(M|M; c) - Q(M|R; c): tr(0) » T(0) (188) 

can be turned into a vector fibration by requiring that the diagram 

K Q(M|R;c) 

TQ(M|R; c) — '—1 _ Q(TM|R; c) 



T Q(M|R;c) 



Q(M 



Q(r M |R; c) 



(189) 



;c) 



l; c) ^^^=^^^= Q(M 

be a vector fibration morphism with the morphism kq(m|r ;c ) defined by 

«Q(M|R:c):TQ(M|R;c) -» Q(TM|R; c): tT(0) ~ q(c,w) (190) 

in terms of a curve 

r: R - Q(M|R; [to, ti]): s h- • ) (191) 

and and a curve 

tu:R->TM:ii->tx(-,t)(0). (192) 

both derived from the same mapping R x R ^ M. Note that kq( . ir ;c ) is a natural transformation 
from the composition TQ( • |R; c) to the composition Q(T • |R; c) of functors T and Q( • |R; c) composed 
in reverse order. 

Covectors are equivalence classes of pairs (q, F) of a configuration q G Q(M|R; c) and a differentiable 
function F on Q(M|R;c). Pairs (q,F) and (q',F') are equivalent if q' = q and 

D(F'or)(0) = D(FoT)(0) (193) 

for each differentiable curve T. The space of cotangent vectors will be denoted by T*Q(M|R; c). The 
equivalence class of (q, F) is denoted by dF(q). The mapping 7Tq( M | R;c ) : T*Q(M|R; c) — > Q(M|R; c) is 
a vector fibration. There is an alternative reprsesentation of covectors. A covector can be represented 
as an equivalence class of triples (c,p, /), where c is a current in R and p and / are curves p: R — > T* M 
and /: R — > T*M such that ttm ° f — P- Triples (c,p, f) and (c',p', /') are equivalent if c' = c 
and 



(f',w) + (p',w) = - / (/,«;> + / <p,w) (194) 
c Joe Jc Joe 

for each mapping w.M. — > TM such that tm ° w — itm ° /• The space of equivalence classes will be 
denoted by P(T*M|R; c) and ^c,p, f) will denote the class of (c,p, /). The mapping 

P(7r M |R; c): P(T*M|R; c) -> Q(M|R; c)S c,p, f) h-» q(c, 7T M ° /) (195) 

is a vector fibration dual to the fibration Q(tm|R; c) with the pairing 

( , }:P(T*M|R;c)x Q(M | R:c)) Q(TM|R;c)^R:((c,p,/),q(c, W ))^- f (f,w) + f (p,w). (196) 

Jc Jdc 
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The diagram 



P(T*M|R;c) 



a Q(M|R;c) 

— — ^ T*Q(M|R; c) 



P(tt m |R; c) 



7I Q(M|R;c) 



(197) 



:Q(M|R;c) 



Q(M|R; c) ^^^=^^^= 
is a vector fibration morphism dual to (189). The mapping Q!q(m|R;c) is defined by 

aQ(M|R;c)( ( c,p, /)) = dF(q(c,7)), 
where the function F and the curves p and / are obtained from the constructions 

F: Q(M|R; c) — ► R: q(c, 7) ^ ^ o t 7 , 

p = 3V t 2 7, 
/ = Sip o t7 

in terms of a mapping </?: TM — ► R and a curve 7: R — > M. 

For a Lagrangian mechanical system the dynamics is contained in the variational principle 



/ (d\,w) = - [ (f,w)+ L M- 
Jc Jc Joe 



(198) 



(199) 

(200) 
(201) 



(202) 



There are two important types of currents in R: an interval c = [to,t\] and a Dirac current c = 
St- If c = [to,ii], then the configuration space Q(M|R; [to, ti]) is the space introduced earlier. If 
c = St then Q(M|R;c) = TM. The equivalence class of (^,7) is the velocity t7(t). Since the 
space Q(M|R; St) is a differential manifold differentiable functions and curves are well defined. Spaces 
TQ(M|R; c) and Q(TM|R; c) are both equal to TTM. Fibrations tq(m|R;c) and Q(tmR; c) assume the 
form ttm:TTM — > TM and Ttm:TTM — > TM respectively. The diagram (189) is in this case the 
familiar isomorphism 



TTM' 



km 



->■ TTM 



TM: 



Tr M 
= TM 



(203) 



of vector fibrations [DA]. The space is the cotangent bundle T*TM. The space T*Q(M|R;c) is the 
fibre product 

TT*M x M T*M = {(p, /) e TT*M x T*M; 7T M (/) - *m(t t * m (p))} . (204) 
The equivalence class of (5 t ,p,f) is the pair (tp(t),f(t). The pairig (196) assumes the form 

( , ):(TT*M x M T*M) x TM TTM ^R:((p,f),w) » (p,w) T - (f,T TM (w)). (205) 
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The domain of this pairing is the set 



{((p, f),w) G (TT*M x M T*M) x TTM; T^ M (p) = Tr M (w) } 
and the pairing ( , ) T is the tangent pairing defined in [DA]. The diagram 



- M 



T7T M 



■ T TM 



71"TM 



(206) 



(207) 



■JM 



is a morphism dual to the vector fibration isomorphism (203). The better known version 



T*TM 



7I"TM 



(208) 



TM: 



TM 



of this morphism is obtained by setting / = and using the tangent pairing 

( , ): TT*Af xjm TTM — > R: (p, w) ^> (p, w) J 

in place of (205). 

The infinitesimal limit of dynamics is formulated in [DA]. 



(209) 



16. Configuration spaces for the statics of continuous media. 

16.1. The configuration of a body, strain. 

Let M be the material space and N the physical space of a medium. The physical space is a 
Euclidean affinc space of dimension 3. By treating this space as a Riemannian differential manifold we 
simplify the presentation of the geometrical framework for statics. The material space is a differential 
manifold. The set of local diffeomorphisms 7: M — > N each defined on an open subset Dom(7) C M 
will be denoted by Q(N\M). Elements of this set are called placement maps. The 1-jet prolongation 
j 1 7: M — > J 1 (A r |M) of a placement map is an important ingredient in the definition of strain. This jet 
prolongation is equivalcntlcntly represented by the tangent map T7. A body is represented by an odd 
current c of dimension m = dim(M) and compact support Sup(c). The vector space of such currents 
will be denoted by CM. We consider commutative diagrams 



J^JVIM)- 



A m T M 



Cl(JV|M) 



(210) 



M ■ 



■M 



where TT m M'- A m T*M — > M is the fibration of odd m-covectors in M and <7i(n\m)'- J 1 (^V|M) — > M is 
the jet-source projection. The space of such diagrams will be denoted by %(N\M). Each diagram is 
well represented by the mapping k: J^iVIM) — > A m T*M. 
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Let CMxQ(7V|M) denote the set 

{(c,7) G CM x Q(JV|M); Sup(c) C Dom( 7 )} . (211) 

In this restricted product set we introduce an equivalence relation. Pairs (c, 7) and (c', 7') are equivalent 
if c' = c and 

/ /cojV = / Koj^ (212) 
Jc Jc 

for each k G %(N\M). The composition k o j 1 7 can be integrated since as a section of the fibration 
7TmM:A m T*M — > M it is an odd m-form on M also called a scalar density. The equivalence class 
of (c, 7) denoted by q(c, 7) is a configuration of the body c. The set of all equivalence classes is 
denoted by Q(N\M). There is the natural projection 7 (tv|M) : Q(N\M) — > CM:q(c, 7 ) c. A fibre 
— T(jv|m) ( c ) ^ s tne configuration space of the body c and the set Q(N\M) is the configuration 

bundle. 

16.2. Vectors and covectors. 

Each element k G 3C(iV|M) induces a function 

n: Q(N\M) — > R: q(c, 7) 1— ► ^ KojV (213) 

Functions constructed in this way are considered differentiable. The space of differcntiablc functions 
on Q(N\M) will be denoted by K(N\M). 

The tangent bundle of the set CM can be identified with the product CM x CM since this set is a 
vector space. 

We consider the set Q(N\M) of differentiable mappings x:txM->JV such that for each s G M 
the mapping x(s, •) is a local diffeomorphism defined on an open subset Dom(x(s, •)) C M. With each 
element (c, 5c, x) of the restricted product 

{{c,5c,x) G CM x CM x Q(N\M); V seR Sup(c + s<5c) C Dom(x(«, ■))} (214) 
we associate the mapping 

g(c, 5c, x): K -> Q(JV|M): s h- q(c + S <5c, x («, ■))■ (215) 

Mappings constructed in this way will be considered differentiable curves. The set of differentiable 
curves will be denoted by G(N\M). 

Tangent vectors in Q(JV|M) are defined as equivalence classes of differentiable curves. Curves g and 
g' are equivalent if 

g'(0)=g(0) (216) 

and 

D( K og')(0) = D(Kog)(0) (217) 

for each differentiable function k. The equivalence class of g will be denoted by tg(0). The space of 
tangent vectors will be denoted by TQ(N\M). It will be shown that the fibres of the projection 

r (JV | M) :TQ(7V|M) - Q(N\M): tg(0) ^ g(0) (218) 

are vector spaces. 

Covectors in Q(N\M) are equivalence classes of pairs (q, k) G Q(N\M) x K(iV|M). Two pairs (q, k) 
and (q', k') are equivalent if q' = q and 

D(K'og(c,Jc,x))(0) = D(#sog(c,<yc,x))(0) (219) 

for each differentiable curve g:R — > Q(N\M) such that g(0) = q. The equivalence class of (q, k) will 
be denoted by d/c(q). The bundle of covectors will be denoted by T*Q(iV|M). Fibres of the fibration 

7T(iv|M): T*Q(N\M) -» Q(JV|M): d/s(q) h-» q (220) 

are vector spaces. 
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16.3. Vertical vectors and covectors. 

A vertical differentiable curve is a mapping 

g(c, x): R - Q(N\M): s ~ q(c, •)) (221) 
associated with an element (c, x) of the restricted product 

{(c, X )eCMxQ(iV|M); V seK Sup(c) C Dom( X ( S , ■))} • (222) 
Vertical curves g and g' are equivalent if 

g'(0) = g(0) (223) 

and 

D( K og')(0) = D(/ t og)(0) (224) 

for each differentiable function k. Equivalence classes are vertical vectors. The equivalence class of a 
vertical curve g will be denoted by tg(0) and the space of vertical vectors will be denoted by VQ(iV|M). 
Fibres of the fibration 

t v {n]m) :VQ(N\M) - Q(AT|M):tg(0) ~ g(0) (225) 

are vector spaces. 

The dual bundle is the set V*Q(iV|M) is the set of equivalence classes of pairs (q, k) £ Q(iV|M) x 
K(N\M). Two pairs (q, k) and (q', k') are equivalent if q' = q and 

D(#s' o g(c, 5c, X ))(0) = D(k o g(c, 5c, x ))(0) (226) 

for each vertical differentiable curve g: M — > Q(N\M) such that g(0) = q. This set is isomorphic to the 
quotient of the bundle T Q(JV|M) of covectors by the polar V°Q(iV|M) of the subbundle of vertical 
vectors. We have the fibration 

^ JV | M) :V*Q(7V|M)^Q(7V|M):tg(0)^g(0) (227) 



16.4. Virtual displacements, force and stress. 

Let 7: M — ► N be a placement map. A commutative diagram 

*■ TiV 




(228) 



is called a deformation of 7. A deformation is well represented by the mapping Sj. Let D{N\M) be 
the space of deformations and let CM~xD(N\M) denote the restricted product set 



{(c,*7) G CM x Q(N\M): Sup(c) C Dom(5 7 )} . 
We consider commutative diagrams 



(229) 



TN ■ 



A m T*M 



TTmM 



(230) 



TV' 



7" 
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and 



TN 



T M 



TN 



Tm-1M 



(231) 



N ■ 



M 



where 7 is a placement map and Tr m -iM- f\ m ~ l T* M — > M is the fibration ol odd (m — l)-covectors 
in M. The diagram (230) represented by the mapping tp is called a force field associated with the 
placement map 7. The diagram (231) is represented by the mapping t and is called a stress field 
associated with the placement map 7. The space of force fields will be denoted by 7{N\M) and the 
space of stress fields will be denoted by 7{N\M). We will use the symbol CMx(7(N\M)x7(N\M)) 
to denote the restricted product composed of triples (c,tp,r) E CM x 7{N\M) x 7(N\M). The force 
field ip and the stress field t in a triple (c, p,r) 6 CMx (5F(7V|M) x7(N\M)) are associated with the 
same placement map 7, and Sup(c) C Dom(7). 

In the set CMxD(N\M) we introduce an equivalence relation. Two pairs (c, £7) and (c',^') are 
equivalent if c' = c and 



tp o $7 



T O £7' = f O Sj + T O Sj 

dc Jc Jdc 



(232) 



for each (c, tp, r) G CM x 3(N\M) x T(iV|M). The equivalence class of (c, ^7) will be called a displace- 
ment of the body c and will be denoted by d(c, £7). The space of all displacements will be denoted by 
D(N\M). There is a dual equivalence relation in the set CMx(J(N\M)x7(N\M)). Triples {c,<p,r) 
and (c', <p',t') are equivalent if c' = c and 



/ tp' o 5-y + / t' o Sj = / tp o 5"f + / t o Sj 
Jc Jdc Jc Jdc 



(233) 



for each (c, £7) e CM xQ(iV|M). The equivalence class of (c, tp, r) will be called a load applied to the 
body c and will be denoted by l(c, tp, t). The space of loads will be denoted by L(N\M). 

16.5. The Euler-Lagrange mapping. 

Fibrations 



^(N\M) x TiV 



T1(N\M)*(tn) 



(234) 



^(N\M) 



and 



J 2 (iV|M) x TAT 

( t 2(JV|M). t n) 



T~2(N\M) {t~n) 



(235) 



J 2 (7V|M) 
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are the inverse images of 



TN 

TN 

N 



(236) 



under the jet-target projections ti(tv|m) : J 1 (iV| A^T) — > AT and T 2 (N\My J 2 (A^|M) — > TV respectively. The 
projections ti(jv|m)*(t/v) and T2(jv|m)*(t/v) are obtained by restricting the projections 



and 

to ^(N\M) 



Wi^N\M)- i\N\M) xTN^ i\N\M) 



pri2 {N \ M )-f{N\M) x TN -> J 2 (Af|M) 



(237) 



(238) 



TAT and J 2 (iV|M) 



TAT respectively. Both inverse images are vector 



( T l(N|M): T lv) (T2(iV|Af) ! r iv) 

fibrations. Fibres {thjm\m * { t n)) ~ 1 {\ 1 l{ x )) and (7"2(jv|M)*(Tiv)) _1 G 2 7( a; )) are both isomorphic to the 
fibre _1 (7(2;)). 

The Euler-Lagrange mapping is a morphism 



J 2 (Ar|M) x TAT 

(T2(N\M) ,Tn) 



T2(N\M)* (tn) 



£dylt 



■ A m T*M 



J 2 (A^|M) 



°2(jV|M) 



(239) 



■ M 



of vector fibrations. The stress-strain mapping is a morphism 

3>d y u 



J^ATIM) x TN 

( t 1(N{M); t n) 



Tl(N\M) (i~n) 



A m - 1 T*M 



1"m-lM 



JH^IM) 



Cl(AT|M) 



(240) 



of vector fibrations. 

A vertical curve in J 1 (A r |M) is a mapping 7:R — > J 1 (A^|M) such that the composition ct\in\m) 7 
is a constant mapping. The image of the composition k o 7 of a vertical curve 7 with a mapping k 
in the commutative diagram (210) is included in the vector space (7r m M) _1 (o'i(jv|M)(7(0)) C R m M. 
The derivative D(k o 7 )(0) is well defined. The vertical differential d V K: VJ^A^jM) i? m M of the 
mapping k is defined in terms of the directional derivative dy K(tj(0)) — D(k o 7)(0). The diagram 



VJ^ATIM) 



dyK 



->- A m T*M 



T P(N\M) 



J^A^IM) 



0"1(JV|M) 



(241) 



->■ M 
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is a morphism of vector fibrations. 

With a function a:M^lwe associate the diagram 

F(a) 



• VJ fe (iV|M) 



VJ fe (iV|M) - 

T J fc (Af|M) r j*(JV|M) 

J fe (iV|M) = ] k (N\M) 

The first row is the mapping 

F(ay.yj k (N\M) ^yj k (N\M):^ k x(0,x ) ^ ^ k X a (0,x ) 

with the mapping x a :IR x M — ► N defined by x a ( s , x ) = x( s ( a ( x ) ~ a ( x o)),x). 
The diagram (242) acts on a diagram 



yj k (N\M) 

T ] k (N\M) 



— A m T*M 



J (N\M) ^— M 



producing the diagram 



lF(a)M 

VJ fe (iV|M) — A m T*M 



'j fc (iV|M) 



where the mapping \F( a )l JL IS defined by 



&1(N\M) 

^(N\M) —— ^ M 



The inverse image 



^(N\M) x T*M 



C1(JV|M)*(7I"m) 



^(N\M) 



diagram is a vector fibration and the mapping 



i F( .)M: (^(N\M) x T*M) x J1{N]M) y^(N\M) -» A m T* M : (da(x) , w) 

("1(K|M) ,7Tm) 

is bilinear. It induces a morphism 
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y^(N\M) 

r J Vl (Af|Af) 



■ JM <g> A m T*M 



&1(N\M) 



(249) 



M 



of vector fibrations. The mapping \p\i in the upper row is characterized by 

(da(x) ® c, ifm(w)} = (c,i F ( )/i(u;)) 



(250) 



with a:M^R,x6M,iD£ TJ 1 (A''|M, a;) C VJ^A^M), and c G A m T x M. We are using in (248) the 
symbol Xji(jv|m) for the fibre product of two vector fibrations. The symbol ir in (249) has an obvious 
meaning. The diagram 



VJ^TVjM) ■ 



l F fJt 



l T*M 



T ^(N\M) 



^(N\M)- 



&1(N\M) 



(251) 



->■ M 



is obtained from the diagram (249) with the use of the isomorphism of TM®A m T*M with A m_1 T*M. 

The mapping ip^i is vertical with respect to the projection t\ (jv|m) m the sense that if Tri(jv|M) ( w ) = 
0, then ipfi(w) = 0. This property makes it possible to construct a mapping 



9n:^(N\M) x TN 

( T l(N|jV/), T iv) 



A r ' 



-l-r* 



T M 



(252) 



characterized by 7'iJ,(j,v) = ip[i(w), where w is a vertical vector such that Tji(jv|m)( w ) = 3 ano - 
~^t\{n\m){w) = v. The diagram 



P(N\M) x TN 

( T l (JV|M) , T iv) 



Tl(jV|M) (T/v) 



7 H 



A m - 1 T*M 



TTm-lM 



J : (^V|M) 



Cl(JV|M) 



(253) 



M 



is a morphism of vector fibrations. 

Let S-y: M — > TiV be a deformation of a placement map 7: M — > A" . The jet prolongations] 1 ^: M — ► 
J 1 (TiV|M) and j 2 <57: M — » J 2 (TiV|M) are transformed in deformations Jj^ = k 1,1 (at| M ) oj 1 ^ and 
S] 2 j = h 2,1 (n\m) j 2 ^7 of the jet prolongations j*7 and j 2 7 respectively. Mappings 7, £7, <5j 1 7, and 
(5j 2 7 can be constructed from a mapping x:Kx M: N by the following definitions 



7: Af — > TV: a; 1 — > x(0,a:), 

^7: M -> TN: x 1 ► tx(-, a)(0), 
flS: Af -> VJ 1 (N\M): x ^ j M x(0, x), 



(254) 

(255) 
(256) 
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and 



(5j 2 7 :M -» VJ 2 (A|M):x j 1,2 x(0, a;), 



(257) 



The total differential is a differential operator applied to linear morphisms from VJ 1 (iV|M) to 
/\m-i-|-*M. The total differential of a morphism v:M] 1 (N\M) -> A m_1 T*M is a linear morphism 
dr^: VJ 1 (A''|M) — > A m_1 T*M characterized by drv o <5j 2 7 = d(zv o <5j 1 7) for each deformation 5-f of a 
placement 7. 

16.6. The internal energy. 

The internal energy density of a material is represented by a mapping u: J 1 (iV|M) ~> A m T*M. The 
internal energy of a configuration q(c, 7) of a body c is the integral 



The internal energy is a function on Q(N\M). 

Variational principles of statics require the definition of the internal energy of a configuration. The 
form of the internal energy varies with the type of the medium. The study of the variational principles 
for different types of media is underway. 
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